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Abstract 

An observation of the A” — » >Jj(2S )A° decay and a comparison of its branching fraction with 
that of the A° —» J/ip A 0 decay has been made with the ATLAS detector in proton-proton 
collisions at yfs = 8TeV at the LHC using an integrated luminosity of 20.6 fb _l . The J/ip 
and ip{2S) mesons are reconstructed in their decays to a muon pair, while the A 0 —» pn~ 
decay is exploited for the A 0 baryon reconstruction. The A? baryons are reconstructed with 
transverse momentum /; T > lOGeV and pseudorapidity |//| <2.1. The measured branch¬ 
ing ratio of the A° —» if/(2S )A° and A^ —» decays is T(A^ —> il/(2S)A°)/T(A° b —> 

J/i[f A 0 ) = 0.501 ± 0.033(stat) ± 0.019(syst), lower than the expectation from the covariant 
quark model. 


© 2015 CERN for the benefit of the ATLAS Collaboration. 

Reproduction of this article or parts of it is allowed as specified in the CC-BY-3.0 license. 




Contents 

1 Introduction 2 

2 The ATLAS detector, data and Monte Carlo simulation samples 2 

3 Event and A® candidate selection 4 

b 

3.1 Charmonium candidate selection 4 

3.2 A 0 and A 0 candidate selection 4 

3.3 A? reconstruction 4 

b 

4 Measurement of the A 0 branching ratio T(A 0 —> <^(2S)A 0 )/T(A 0 -4 J/ijj A 0 ) 8 

b b b 

5 Summary 10 


1 Introduction 

The A® baryon properties have been extensively studied at the Large Hadron Collider (LHC) [1-7]. The 
decay channel A ( ; J —> J/ij/(jj. + iu~)A°(pn~) 1 has been primarily used by the LHC experiments in these 
studies, although a number of other A° b decay channels have been exploited by the LHCb experiment. 
In particular, a measurement of the differential branching fraction and angular analysis of the rare decay 
A® —> yu + /CA° was performed by LHCb [8,9] following observation of this decay by the CDF experi¬ 
ment [10] at the Tevatron collider. However, no results for the decay mode A? —» iIj(2S )A° have yet been 
reported, although a measurement of the decay properties would be useful for verification of theoretical 
predictions [11]. 

The A° —» (Jj(2S )A° branching fraction should be of the same order as that of the decay A® —> J/if/A° as 
suggested by the branching fraction values of the B°, B + and B° s meson decays to >Jj(2S) / J / ip and either 
a pseudoscalar ( K °, K + , if) or vector (K*°, K* + , <p) meson. The branching ratios of such B meson decays 
to iJ/(2S)X and J/ipX are within the 0.5-0.8 range [12], and are generally reproduced by factorisation 
calculations [13]. The only available theoretical calculation of the branching ratio of the A^ —> iJj(2S )A° 
and A*’ —> J/tf A 0 decays, performed in the framework of the covariant quark model [14], predicts 0.8 
with an uncertainty of approximately 0.1 [11]. 

An observation of the A® —> if/(2S)A° decay and a measurement of the branching ratio of the A® —> 
if(2S)A° and A^ -4 J/ifA° decays is reported in this Letter. The J/ip and <//( 2S) mesons are reconstructed 
in their decays to a muon pair, while the A 0 —> pn decay is exploited for the A 0 baryon reconstruction. 
The A» baryons are reconstructed with transverse momentum p-\ > lOGeV and pseudorapidity \tj\ < 
2 . 1 . 


2 The ATLAS detector, data and Monte Carlo simulation samples 

A detailed description of the ATLAS detector can be found elsewhere [15]. A brief outline of the com¬ 
ponents most relevant to this analysis is given below. 


1 Hereafter, charge conjugation is implied, unless explicitly stated otherwise. 
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The ATLAS inner detector (ID) has full coverage- in (f>, covers the pseudorapidity range |//| < 2.5 and 
operates inside an axial magnetic field of 2T. It consists of a silicon pixel detector (Pixel), a silicon 
microstrip detector (semiconductor tracker, SCT) and a transition radiation tracker (TRT). The inner- 
detector barrel (end-cap) parts consist of 3 (2 x 3) Pixel layers, 4 (2 x 9) double-layers of single-sided 
SCT strips and 73 (2 x 160) layers of TRT straws. The ATLAS muon spectrometer (MS) covers the 
pseudorapidity range |//| < 2.7. It consists of precision tracking chambers, fast trigger detectors and a 
large toroidal magnet system generating an average field of 0.5 T in the barrel region (|//| < 1.05) and 1 T 
in the end-cap regions (1.05 < H < 2-7). 

The ATLAS detector has a three-level trigger system [16]: the hardware-based Level-1 system and the 
two-stage High Level Trigger (HLT). For this measurement, dimuon triggers are used. At Level-1, the 
dimuon triggers search for patterns of MS hits corresponding to dimuons passing various pj thresholds. 
Since the rate from the low-px dimuon triggers was too high, prescale factors were applied to reduce their 
output rates. The data sample used in this analysis was collected using three dimuon triggers with pj 
thresholds of 4 GeV for both muons, 4 GeV and 6 GeV for the two muons, and 6 GeV for both muons. At 
the HLT, the dimuon triggers used require muons with opposite charges and dimuon mass in the range 
2.5 < m(p + p~) < 4.3 GeV. 

This analysis uses 20.6 fb' 1 of proton-proton collision data with a centre-of-mass energy of 8 TeV recor¬ 
ded by the ATLAS detector at the LHC in 2012. The uncertainty on the integrated luminosity is ±2.8%. It 
is derived following the same methodology as that detailed in [17]. The event sample is processed using 
the standard offline ATLAS detector calibration and event reconstruction code. There are typically a few 
primary vertex candidates in each event due to multiple collisions per bunch crossing. Only events with at 
least four reconstructed tracks with px > 0.4 GeV and at least one reconstructed primary vertex candidate 
are kept for further analysis. 

To model inelastic pp events containing A <3 —> J/i//(p + p~) A 0 , A® —» i//(p + p~) A 0 , B° —» J/i//{p + p~)K ° 
or B° —> iJ/(ju + p-)Kf decays, 2 3 four large samples of Monte Carlo (MC) simulated events are prepared 
using the PYTHIA 8.1 [18] MC generator. The B° MC samples are needed to control reflections from 
B° decays to the A° signal distributions. The generation is based on leading-order matrix elements for all 
2 —> 2 QCD processes. Initial- and final-state parton showering is used to simulate higher-order processes. 
Generated events with both muons from 7/t/r or t//( 2S ) decays having transverse momenta above 3.5 GeV 
and pseudorapidities within ±2.5, and, for A° MC samples, with the A 0 transverse momentum above 
1 GeV are passed through a full simulation of the detector using the ATLAS simulation framework [19] 
based on GEANT4 [20,21] and processed with the same reconstruction program as used for the data. An 
emulation of the three triggers used for the data collection is applied to the MC samples. The angular 
decay distributions of the Ajj —> J/i//(p + p~)A°(pn~) decay are modelled using the helicity amplitudes 
measured by ATLAS [2], For the A® —> ip(p + p~)A {] (pn~) decay, the helicity amplitudes are set to the 
predicted values [11]. 


2 The ATLAS coordinate system is a Cartesian right-handed system, with the coordinate origin at the nominal interaction 
point. The anti-clockwise beam direction defines the positive z-axis, with the x-axis pointing to the centre of the LHC ring. 
Polar ( 6 ) and azimuthal (0) angles are measured with respect to this reference system. The pseudorapidity is defined as 
rj = - lntan(0/2). 

3 In this Letter, i//(2S) is referred to as ip when its decay channel is indicated. 
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3 Event and A 0 candidate selection 

b 

3.1 Charmonium candidate selection 

Events are required to contain at least two muons identified by the MS with tracks reconstructed in the ID. 
The reconstructed muons are required to match the muon candidates identified by the trigger. The muon 
track parameters are taken from the ID measurement alone, since the MS does not significantly improve 
the precision in the momentum range relevant for the charmonium measurements presented here. To 
ensure accurate measurements, each muon track must contain at least six SCT hits and at least one Pixel 
hit. Muon candidates satisfying these criteria are required to have opposite charges and a successful fit to 
a common vertex with A 2 /AW < 10, where A 2 is the fit quality with the number of degrees of freedom 
/Vjof = 1. Events with m(ju + p ~) values within ±200 MeV intervals around the //</r and )Jj( 2 .S ') world 
average masses [12] are used to search for A 0 —> pjx~ candidates. 


3.2 A 0 and A 0 candidate selection 

In all events with J/iJ/ or if/(2S) candidates, pairs of tracks from particles with opposite charge are com¬ 
bined to form A 0 candidates. Each track is required to have at least one Pixel or SCT hit. Only pairs 
successfully fitted to a common vertex with W/AW < 5 are kept. The track with larger momentum 
is assigned the proton mass hypothesis since the proton always has a larger momentum than the pion 
for A 0 baryons with momenta larger than 0.3 GeV. To suppress combinatorial background the following 
requirements are used: 

• Pt(p) >1-7 GeV. 

• |zoO)l < 25 mm, where zo(p) is the proton longitudinal impact parameter with respect to the dimuon 
vertex. MC studies show the requirement produces no loss of signal. 

. L^(A 0 ) > 7 mm, where L^(A°) is the transverse decay length 4 of the A 0 candidate measured 
from the beam line. 

Events with m{pn~) values within an interval of ±20MeV around the A 0 world average mass [12] are 
kept for further analysis. 


3.3 A 0 reconstruction 

b 

Tracks of the selected charmonium and A 0 candidates are simultaneously refitted with the dimuon and di¬ 
hadron masses constrained to the world average masses of //i/f ( mj y^) or i//(2S ) (rn^S)) and A 0 (m A o) [12], 
respectively. The combined momentum of the refitted A 0 track pair is required to point to the dimuon 
vertex. To control B° reflections to the A® signal distributions, a B° decay topology fit is also attempted 
for each track quadruplet successfully fitted to the A 4 topology, i.e. the pion mass is assigned to both 
hadron tracks and the dihadron mass is constrained to the world average mass of K® [12]. To suppress 
combinatorial and B° backgrounds the following requirements are used: 


4 The transverse decay length of a particle is the transverse distance between the primary or production vertex and the particle 
decay vertex projected along the transverse momentum of the particle. 
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Figure 1: The invariant mass distributions m(J/tf/A°) (left plot) and m(J/i//A °) (right plot) for selected A° b and A { b 
candidates, respectively. The solid histograms represent the fit results (see text). The A'J signals (dashed lines) and 
the B° reflections are also shown. 




m(JAi»A°(A )) [MeV] 


m(J/y K°) [MeV] 


Figure 2: The invariant mass distributions for the combined sample of the selected A ( b and A® candidates obtained 
after their fits to the A b —> J/if/ A 0 (left plot) and B° —> J/if/K ° (right plot) topologies. The solid histograms 
represent fit results (see text). The A ° b and B° signals and their mutual reflections are also shown. 
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Figure 3: The invariant mass distributions for the combined sample of the selected and A° candidates obtained 
after their fits to the A° —> if/(2S)A° (left plot) and B° —> tf/(2S)K° (right plot) topologies. The solid histograms 
represent fit results (see text). The A y and B° signals and their mutual reflections are also shown. 


• ^(A^/ZVjof < 3, where y 1 is the quality of the fit to the A? topology with Vdof = 6. 

• L xy ( A 0 ) > 10 mm, where L xy ( A 0 ) is the transverse decay length of the refitted A 0 vertex measured 
from the A? (dimuon) vertex. 

• p T (A°) > 2.5 GeV. 

• pj{n~) > 0-45 GeV. 

• r(A°) > 0.35 ps, where r(A°) = L xy { A°) • m A o/p^A^) is the A° proper decay time, L xy (A° b ) is the 
transverse decay length of the A® vertex measured from the primary vertex and m A o is the A® world 
average mass [12]. The primary vertex candidate with at least three tracks and the smallest value 
of the three-dimensional impact parameter of the A° candidate is selected as the actual primary 
vertex. 

• P{ A°) > f»(fi°), where P(A°) and P(B°) are the x 2 probabilities of the quadruplet fits with A® and 
B° topologies, respectively. 

The muon transverse momenta and pseudorapidities are required to be in the ranges with high values of 
the trigger and reconstruction acceptances: 


Pt(m ± )> 4 GeV, 1^)1 <2.3. 
The kinematic range of the A b measurement is fixed to 


p T (A° b ) > 10 GeV, |i/(A°)|<2.1. 

The invariant mass distribution m(J/ij/ A 0 ), calculated using track parameters from the A b topology tits, 
is shown in Figure 1 separately for the selected A° and AJJ candidates. Clear signals with similar size 
are seen in the two distributions around the world average mass of the A® baryon. Figures 2 and 3 
show the m(J/i// A 0 ) and m(if/(2S )A°) distributions for the combined sample of the A ( b and A® candidates. 
The invariant mass distributions and m(iJj(2S )K^) from the B° topology tits are also shown. 
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Clear signals are seen in the m(J/if/ A 0 ) and m(ip(2S)A°) distributions 5 around the world average mass 
of the A® baryon. There arc also signals in the m(J/ifsK®) and m(iJ/(2S )K®) distributions near the world 
average mass of the B° meson [12]. The B° signals are smaller than the signals due to the selection 
requirements. 

The A 0 ) and m(J/ipK^) distributions are simultaneously fitted to sums of signal and two-component 
background distributions. The signals are described by modified Gaussian functions [22]. The modified 
Gaussian function is defined as 


Gauss mod oc exp[-0.5 ■ * 1+1 /(i+o. 5 .*) L 

where x = |(m - mo)/cr|. This functional form, introduced to take into account the non-Gaussian tails 
of resonant signals, describes both data and MC signals well. The signal position, mo, and width, <x, as 
well as the number of the signal events arc free parameters of the fit. The non-resonant backgrounds in 
the distributions are described by independent exponential functions. The mutual B° and A° reflections 
are described by MC templates normalised to the numbers of B {) and A'J hadrons obtained in the fit. The 
reflection normalisations are corrected for small losses (2-6%) of A d and B° hadrons that passed the 
A® reconstruction but failed the B° reconstruction. The corrections are obtained using MC simulation. A 
similar fit is performed for the m(ip(2S )A°) and m( t//( 2S )K ®) distributions. In the analysis of the combined 
A® and A® samples, the ratio of the MC A*’ and A” events is set to the data ratio obtained in the separate 
A'J —> J/ipA° and A® —> J/ipA° fits (Figure 1). The A° and A° fitted yields are 3523 ± 89 and 3414 ± 92, 
respectively, providing the ratio 1.03 ± 0.04(stat). 

The results of the fits for the combined A® and A (5 samples are summarised in Table 1. The Ajr mass 
values obtained from the fits of the m(J/ip A 0 ) and m(ip(2S )A°) distributions agree with each other and 
with the world average A° mass value [12]. The signal widths are different, reflecting the difference in 
charmonium masses in the two decay channels, in agreement with the MC expectations. The quality, 
yA/Adof, of the A® —> J/ifj A° and A <5 —> iJ/(2S)A () signal fits are 1.0 and 1.1, respectively. 



A° b -» J/if/A® 

B° -> J/ifjKl 

A° -» tp(2S) A 0 

B° -> ip(2S)K° s 

j\J ■ 
iV sig 

6940 ± 130 

854 + 84 

603 ± 38 

124 ± 28 

m sig [MeV] 

5620.4 ± 0.4 

5274.7 ± 2.3 

5618.2+ 1.2 

5272.4 ± 4.9 

cr sig [MeV] 

19.7 + 0.5 

19.2 + 2.2 

14.3 ± 1.1 

16.7 + 4.1 


Table 1: The numbers of signal events, A s j g , signal masses, m sl g, and signal widths, cr s i g , obtained by the fits (see 
text). Only statistical uncertainties are shown. 


To verify that the observed A*’ signals correspond to the A° —> J/ipA° and A*,’ —> iJj(2S )A° decays the 
signal reconstruction is repeated with only one mass constraint for either the dimuon or the dihadron track 
pair in the cascade fit and the A*’ mass is calculated using the mass-difference method. In the case that 
the dihadron mass is fixed to the A° mass, the A® mass is calculated as m(ji + n~ A 0 ) - m(p + p ~) + m.//,/, 
(m{p + p~ A 0 ) - m(jj. + n ~) + m^s)) for m(ju + p~ ) < 3.4 GeV (m(p + p~) > 3.4 GeV). When the dimuon mass 
is fixed to the J/ip (iJ/(2S j) mass, the A? mass is calculated as m(J/i/jpn~) - ni(pn~) + m A n (m(ip(2S)pn^)~ 
m(pK~)+m A o). In both cases clean A 15 signals are reconstructed with numbers of signal events compatible 

5 Studies with MC simulated events show that the fraction of reconstructed —> J/if/A 0 decays which can contribute to the 

reconstructed A® —> i//(2S )A° signal is ~ 10 5 . 
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Figure 4: The m(jj. + p~) distributions for A° b —» J/i/zA 0 candidates (left plot) and A° b —> i//(2S)A° candidates (right 
plot) after full selection, without a mass constraint on the charmonium mass in the cascade fit. The spectra are fitted 
with a sum of an exponential function and a modified Gaussian function. 


with those in Table 1. Figure 4 shows the m(p + p~ ) distributions for the A° —> J/ijj A 0 and A® —> t/r(25)A° 
candidates reconstructed with the mass constraint for the dihadron pair and selected within ±3<x s i g around 
the world average A‘| mass. Clear signals from T/i/r and ip(2S) are seen. The m(ju + p~) distributions are 
fitted to a sum of an exponential function describing the background and a modified Gaussian function 
describing the signal. The signal yields are found to be Nj/# = 9770 ± 120 and N^(2S) = 724 ± 45. 
Figure 5 shows the m(pn~) distributions for the A° candidates reconstructed with the mass constraint for 
the dimuon pair and selected within ±3o" s j g around the world average A® mass. Clear signals from A 0 are 
seen. The m(pn~) distributions are fitted to a sum of a threshold function describing the background and 
a modified Gaussian function describing the signal. The threshold function has the form 

d 'y 

A ■ (m — m p — m„- ) ■ exp[C • (m - m p - m n ) + D ■ (in - m p - m„ ) ], 

where m p and m n - are the proton and pion masses, respectively, and A, B, C and D are free parameters. 
The A 0 signal yields are found to be 7710 ± 120 and 702 ± 38 for the A° —> J/if/A° and A® —> >//(2S )A° 
candidates, respectively. The numbers of signal charmonium and A 0 events are larger than the numbers 
of the corresponding A® signal events because the backgrounds are partly due to genuine charmonium 
and A 0 states. 

4 Measurement of the .V' branching ratio 

r(A“ -> tfr(2S)A°)/r(A“ -> JN> a 0 ) 

The numbers of A° signal events in the two decay modes, reported in Table 1, are corrected for detector 
effects and selection efficiencies as N cor = N s i g /Sl, where N COI is the corrected number and J\ is the MC 
acceptance. The MC events with the i//(2S)/J/ifi muons having transverse momenta above 3.5 GeV and 
pseudorapidities within ±2.5, and A 0 transverse momentum above 1 GeV, passed through the detector 
simulation and event reconstruction, are used to correct the numbers of signal events in the fiducial range, 
defined as follows: 














Figure 5: The m(pn~) distributions for A° b —» 7/i/rA 0 candidates (left plot) and —» i//(2S)A° candidates (right 

plot) after full selection, without a mass constraint on the A 0 mass in the cascade fit. The spectra are fitted with a 
sum of a threshold function and a modified Gaussian function. 


Pt(A^) > lOGeV, | 77 (AO)| < 2.1, 

PiUf) > 4GeV, I 17 GOI < 2.3, 
p T (A°) > 2.5 GeV. 

The acceptances are calculated as the ratio of the number of reconstructed A° signal events passing all 
selection requirements in the above fiducial range to the number of A? baryons in the same decay mode 
and fiducial range at the MC generator level. These acceptances are 4.16 ± 0.02(stat)% and 4.30 ± 
0.03(stat)% for the A ( ; J —> J/iJ/A 0 and A® —> iJ/(2S )A° decays, respectively. In the fiducial range, the ratio 
of the corrected numbers of A® signal events in the two decay modes is 0.0841 ± 0.0055(stat). 

Then the numbers are corrected, using generator-level MC samples with no requirements on the and 
A 0 selection, from the above fiducial range to the kinematic range of the A® measurement 


Pt(A^) > lOGeV, |J7(A°)|<2.1. 

The acceptances of the latter corrections are 7.57±0.06(stat)% and 9.61±0.07(stat)% for the A*J —> J/ij/A 0 
and A° —> i//(2.S')A° decays, respectively. Finally, the branching ratio of the two A ( ; J decays is calculated 
as 


T(A° -> ft(2S )A°) _ N COI (A° b -> if/(p, + ^ -> Tt) 

T(A« -» J/if/A°) ~ 7V cor (A° -» //(A(yU + jU~)A°) ' &M2S) -» £+{-)' 

where & is the branching fraction of the corresponding charmonium decay to a lepton pair. In the case of 
7/i/r, the branching fraction S(//i/r —> p + pT) = 0.05961 ± 0.00033 [12] is used. For S(i//(2S) —» l + £~), 
the branching fraction S(iJ/(2S) —> e + e~) = 0.00789 ± 0.00017 is used, assuming lepton universality, 
because it is measured with better precision than in the muon channel, S(i/r(2>S) —> = 0.0079 ± 

0.0009 [12]. 

Five groups of systematic uncertainty sources are considered. The effect of each group on the measured 
ratio, obtained by adding in quadrature the effects of independent sources, is shown in parentheses: 
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• Dependence on the production model (±0.1%). The uncertainty is obtained by 

- varying the MC pj(A®) and |t/(A?)| distributions while preserving agreement with the data 
distributions, 

- varying the MC ratio of A? and A*’ baryons in the range allowed by the separate data fits 
(Section 3), 

- varying the lifetimes of the A 0 and A ^ baryons in the ranges of their uncertainties [12]. 

• Dependence on the A*’ polarisation model (±1.1%). The uncertainty is obtained by varying the 
MC A® —> J/i//(p + p~)A°( pn~) helicity amplitudes in the range of their uncertainties [2], and by 
changing the MC A? —> d/(u + p~)A 0 (p7:~) helicity amplitudes to those measured by ATLAS for the 
A° -» // i//(p + p~)A°(pn ~) decay [2], 

• The uncertainty of the signal extraction procedures (±2.8%). The uncertainty is determined by 
changing the background parameterisations to second order polynomials and by reducing the ranges 
used for the signal tits by 20MeV from either left or right side, independently for the two A {) h 
signals. In addition, the corrections of the reflection normalisations, obtained from MC simulation, 
are varied by half of their values. This uncertainty is affected by statistical fluctuations. 

• The uncertainty originating from the MC statistical uncertainty (±1.3%). 

• The uncertainty of the charmonium branching fractions B(J/ip —> p + p~) and B(4)(2S) —> e + e~) 
(± 2 . 2 %). 

The measured branching ratio of the two A° decays is 
T(A? -» <A(25)A°) 

- 'A -= 0.501 ± 0.033(stat) ± 0.016(syst) ± 0.011(S), 

T(A® —» J/i/zA 0 ) 

where the contributions from the first four groups of systematic uncertainty are added in quadrature. The 
uncertainty due to the uncertainties of the charmonium branching fractions B is quoted separately. The 
luminosity uncertainty, uncertainties of the muon and hadron track reconstruction and the vertexing un¬ 
certainties cancel out in the ratio. The bias in the measured ratio due to contributions from the rare decay 
A? —> p + p~ A 0 is estimated using the LHCb measurement [9] of the rare decay’s differential branching 
fraction to be below 0.5% and thus neglected. Consistent ratio values are found when calculated in bins 
of Pt{A®) or separately for A° and A® baryons. 

The measured ratio lies in the range 0.5-0.8 found for the branching ratios of analogous B meson de¬ 
cays [12]. The only available calculation for the branching ratio of the two A^ decays (0.8 ±0.1 [11]) 
exceeds the measured value. 


5 Summary 

The A° —» iJ/(2S )A° decay has been observed with the ATLAS detector in pp collisions at \fs = 8TeV 
at the LHC using an integrated luminosity of 20.6 fb _1 . The branching ratio of the A^ —> iJ/(2S )A° 
and A® —» J/iJ/A 0 decays has been measured to be T(A® —> i/r(25')A°)/T(A^ —» J/iJ/A () ) = 0.501 ± 
0.033(stat) ± 0.016(syst) ± 0.01 1(B). The ratio falls into the range 0.5-0.8, as found for the branching 
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ratios of analogous B meson decays [12]. The only available theoretical expectation for the branching 
ratio of the two decays (0.8 ±0.1 [11]) exceeds the measured value. 
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